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ABSTRACT 

In this paper we compute the gravitational signal emitted when a white dwarf 
moves around a black hole on a closed or open orbit using the affine model 
approach. We compare the orbital and the tidal contributions to the signal, 
assuming that the star moves in a safe region where, although very close to 
the black hole, the strength of the tidal interaction is insufficient to provoque 
the stellar disruption. We show that for all considered orbits the tidal signal 
presents sharp peaks corresponding to the excitation of the star non radial 
oscillation modes, the amplitude of which depends on how deep the star pen- 
etrates the black hole tidal radius and on the type of orbit. Further structure 
is added to the emitted signal by the coupling between the orbital and the 
tidal motion. 

Key words: gravitational waves — stars: oscillations — relativity — stars:white 
dwarfs 



1 INTRODUCTION 

Recent astronomical observations have shown evidence of interactive processes in action 
between stars and black holes. The ESO Very Large Telescope (VLT) has observed stars 
that move on very elliptic orbits around the supermassive black hole at the center of our 
Galaxy (Schodel et al. 2002), and X-ray satellite CHANDRA has monitored the event of a 
star disrupted and swollen by the giant black hole at the center of the galaxy RXJ1242-11 
(Komossa et al. 2004). These and similar phenomena that take place around black holes 
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generate gravitational waves essentially because of two mechanisms: one related to the time 
variation of the quadrupole moment of the system star-black hole due to the orbital motion, 
the second to the quadrupole moment of the star which changes in time because the star 
is deformed by the tidal interaction with the black hole. The orbital contribution has been 
computed using the Post-Newtonian formalisms (see for instance (Blanchet 2002) and ref- 
erences therein). In this paper we shall focus on the second mechanism: we shall consider a 
white dwarf (WD) moving on assigned orbits around a black hole (BH), we shall compute 
how the stellar shape and structure change along its motion and the gravitational wave 
(GW) signal emitted due to this time-varying deformation. 

The problem of tidal deformation has been widely investigated both in newtonian gravity 
and in general relativity using different approaches. 

In 1977 Press and Teukolsky studied the formation of binary systems by two-body tidal 
capture using a perturbative approach and provided estimates of the amount of orbital 
energy absorbed in the encounter (Press & Teukolsky 1977). 

In the early eighthies a formalism was developed (Carter & Luminet 1982; Carter & 
Luminet 1983; Carter & Luminet 1985) that allows to compute the tidal deformation by 
integrating a set of equations which describe the motion of an element of stellar fluid due 
to the effect of the tidal tensor of the other massive body on the assumption that while 
deforming the star maintains an ellipsoidal form (affine model). This approach, initially 
developed in the framework of Newtonian gravity, was subsequently generalized to general 
relativity when one of the interacting bodies is a Schwarzschild or a Kerr black hole (Marck 
1983; Luminet & Marck 1985). The affine model has been used in several papers addressing 
some relevant astrophysical problems: in (Luminet & Carter 1986; Luminet & Pichon 1989; 
Bicknell & Gingold 1983) the authors investigated whether, in the interaction with a mas- 
sive companion, a main sequence (MS) star can be compressed to such an extent that the 
temperature increase in the core can ignite a detonation of the nuclear fuel (a process which 
may occur in active galactic nuclei having a central massive black hole); the tidal capture 
of a main sequence star and its disruption have been investigated in (Carter 1992) and in 
(Kochanek 1992), respectively; tidal processes occurring in neutron star- neutron star coa- 
lescing binaries before contact has been investigated in (Kochanek 1992), where the phase 
shift induced on the emitted gravitational signal has also been computed; the chemical and 
mechanical behaviour of a WD passing inside the tidal radius of a black hole has been 
studied in (Luminet & Pichon 1989); furthermore, the affine model has also been used in 
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(Kosovichev & Novikov 1992; Diener et al. 1995) to study the tidal capture and interaction 
of a star and a massive BH. 

Tidal effects in close compact binaries have also been studied using different formalisms; 
for instance, by integrating the hydrodynamical equations in 3D, in a form adapted to 
express the evolution of the star principal axes and of other relevant quantities, and including 
radiation reaction and viscous dissipation (Lai, Rasio & Shapiro 1993; Lai, Rasio & Shapiro 
1994; Lai, Rasio & Shapiro 1994; Rasio & Shapiro 1995; Lai, Rasio & Shapiro 1995). An 
alternative approach to the integration of the 3D-hydrodynamical equations has been used 
in (Frolov et al. 1984) to study disruptive and non disruptive encounters between white 
dwarfs and black holes, focusing on the determination of the periastron distance, time delay, 
relativistic precession and tidal forces. 

Further studies on tidal interactions in close binaries are in (Evans & Kochanek 1989; 
Novikov, Pethick & Polnarev 1992; Laguna et al. 1993; Khokhlov, Novikov & Pethick 1993; 
Khokhlov, Novikov & Pethick 1993; Ogawaguchi & Kojima 1996; Marck, Lioure & Bonazzola 
1996; Diener et al. 1997; Fryer et al. 1999; Ferrari, D'Andrea & Berti 2000; Shibata & Uryu 
2001; Gualtieri et al. 2001; Berti k Ferrari 2001; Ivanov & Novikov 2001; Pons et al. 2002; 
Ivanov, Chernyakova & Novikov 2003; Gomboc & Cadez 2005). 

In this paper we use the Carter-Luminet approach to study the tidal interaction of a 
white dwarf and a black hole, with the purpose of computing the gravitational signal emitted 
in this process. We shall integrate the Carter-Luminet equations for open and closed orbits 
choosing a white dwarf mass M* = 1 M and a black hole mass M B h = 10 M . The WD will 
be modeled using a polytropic equation of state (EOS) with adiabatic index 7 = 5/3. The 
gravitational radiation will be computed by using the quadrupole approach and the orbital 
and the tidal contributions will be compared. We shall further discuss how the results change 
for higher values of the black hole mass. 

Our study is based on several simplifying assumptions: 

- The WD mass is assumed to be much smaller than the black hole mass. 

- The black hole does not rotate and the equations of motion of the star, its internal structure 
and the black hole tidal tensor are computed using the equations of newtonian gravity; thus, 
we call "black hole" an object which is basically a point mass star. We plan to extend our 
calculations to general relativity and to rotating black holes in the near future. 

- We neglect tidal effects on the orbital motion. 
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- We neglect viscous effects, an hypothesis which appears to be justified in the case of 
WD-BH binaries (Wiggins & Lai 2001). 

The plan of the paper is the following. In section 2 we summarize the equations we use 
to compute the tidal deformation of the star in the affine model approach, and those needed 
to compute the emitted radiation. In section 3 we discuss the results of our study: for the 
1 M Q — 10 M Q WD-BH binary we compare the orbital and the tidal contributions to the 
gravitational signal, for circular, elliptic and parabolic orbits; we further discuss how these 
results depend on the black hole mass. Conclusions are drawn in section 4. 



2 THE RELEVANT EQUATIONS 

As mentioned in the introduction, the equations we use to describe the tidal deformation of 
the star are fully Newtonian, and have been developed in (Carter & Luminet 1982; Carter 
& Luminet 1983; Carter & Luminet 1985). 

Let us consider a spherical star with mass M* and radius i?* in equilibrium, and be f 
the Cartesian position vector of the generic fluid element with respect to a frame relative to 
the center of mass of the star in this unperturbed configuration. To hereafter, the 'hat' will 
indicate quantities computed in this frame for the unperturbed star. 

We shall assume that the star is composed of a perfect fluid with a polytropic EOS 

P = Kp\ (1) 

M* will move on assigned orbits around a non rotating, massive body M BH ^> M*, to which 
we associate a reference frame Obh- Be X(t) the vector which identifies the position of the 
center of mass of the star with respect to this frame. 

While moving around Mbh, the star is deformed by the tidal interaction, and we shall 
indicate as O* the reference frame relative to the center of mass of the moving star in a 
parallel propagated frame; in this frame, the position of each fluid element will be given by 
the vector f(t). Thus, the position of the generic fluid element of the star with respect to 
Obh will be 

Xi = Xi + n, (2) 
and f will be a function of time, related to f by 

n = qijfj (i, j = 1, 2, 3) (3) 
where qij is a 3 x 3 uniform matrix, i.e. it depends on time but it is independent of r, which 
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implies that the star can assume only ellipsoidal configurations (affine model), is the 
matrix we need to compute to know how the stellar shape changes along the orbit. 

Under the assumption of polytropic EOS, the satisfy the following equations (Carter 
& Luminet 1983) 

^ 2 Jo duK —^^ , (4) 



4ij — CikQkj + 



IL 



where 



Cij 



dXidXj 



is the black hole tidal tensor computed from the Newtonian potential 
GM BH 



E 



X 



In eqs. (4) II is the volume integral of the local pressure 
n = J P-dV, 



and, since 



P = p- \\q\ 



(5) 



(6) 



(7) 



(8) 



where p is the density in the spherical unperturbed configuration and ||g|| = det q, it can 
be written as 



n = K p^dV = K\\q 



1 1-7 J p7~l 



dM 



\Q\ 



1 1-7 



Pdv = n* \ \q\ 



11-7 



(9) 



II* is the value of II in the unperturbed configuration. Furthermore, M.* is the scalar 
quadrupole moment of the initial configuration 

M * = lJ f ^ dM - ( 10 ) 
The function A in eq. (4) is 

A = y/\\S + uI\\, (11) 
where I is the identity matrix, 

Sij = qikQjk, (12) 

and u is the integration variable of the integral appearing in (4). 

If the star is initially non rotating (as we shall assume), or if its initial spin is orthogonal 
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to the orbital plane, the non vanishing components of reduce to (Carter & Luminet 1983) 
( <?n Q12 \ 

(13) 



Qij(t) = 



521 <?22 

o o g 33 y 



2.1 Gravitational radiation 

The gravitational signal emitted by the system will be computed by using the quadrupole 
formalism, according to which 



h T J = n = 0, 3 



where mf T is the source quadrupole moment projected on the TT-gauge 



rrii 



'PijlmJR'ir. 



(14) 



(15) 



and Vijim is the projector tensor. For the system under consideration, the quadrupole mo- 
ment can be evaluated as follows 



m ik (t) = -J / T 00 (t, X j) XiX k d 3 x = f p(t, rj) {X, + n) {X k + r k ) d 3 r, 



(16) 



i.e. 



m ik (t) = J v p(t : r j )X i (t)X k (t)d 3 r + (t, rj ) nr k d 3 r 

+ X i {t)J v p{t ) r j )r k d z r + X k {t)J v p{t,r j )r l d i r ) (17) 

where we remind that the Xi(t) are the coordinates of the center of mass of the star with 
respect to Obh, whereas r k (t) are the coordinates of the element of stellar fluid with respect 
to O*. The first integral is the standard quadrupole moment associated to the orbital motion 

m°l b (t) =X l (t)X k (t) f p(t,r :j )d :i r = M*X l (t)X k (t). (18) 

The second integral is the contribution due to the tidal interaction which, using eqs. (3) can 
be written as 

™>ik f (t) = J v p(t,rj)rir k d 3 r = qu(t)q km (t) J p (t,fj) fif m d 3 f. (19) 
In the last equality we have used eq. (8) and the property that the volume scales as 
d 3 r = \\q\\d 3 r. (20) 
Using eqs. (10) and (12), eq. (19) becomes 

mt f {t) = qu(t)q km (t) [ hf m dM = M*S ik {t). (21) 
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Gravitational waves from tidal interactions 7 

Table 1. The limiting value of /3 above which the tidal interaction between the white dwarf and the black hole becomes 
disruptive is given for the orbits considered in this paper. 



Orbit 


e 


Pcrit 


Circular 





0.78 


Parabolic 


1 


1.05 


Elliptic 


0.25 


0.78 


0.75 0.97 




0.95 


1.03 



The last two integrals in (16) vanish due to symmetry; indeed 

Xi(t) J p (t, Tj) r k d 3 r = X t (t)q km (t) p (t, f 3 ) f m d 3 f = (22) 

since for every fluid element that, in the unperturbed configuration, has coordinates rj = 
(xj, i/j, Zj) there exists one with coordinates (—Xj, —i/j, —Zj). Thus, in conclusion, the 
quadrupole moment of the system is composed of two contributions, the orbital one, m," 6 ( 
eq. 18 ), and that due to the tidal interation, (eq. 21). 



3 RESULTS 

Eqs. (4) have been integrated for a system composed of a white dwarf described by a poly- 
tropic EOS with 7 = | and a more massive companion of mass M B h = 10 M Q . The WD 
central density of the unperturbed configuration is such that M* = 1 M & and R* = 6955 km. 

We have considered circular, elliptic and parabolic orbits, and using eqs. (18 ), (21) and 
(14) we have computed the emitted gravitational signal. 

It is known that if the periastron of the orbit lays within a certain critical radial distance 
from Mbh the star can be destroyed by the tidal interaction. This radius, Rr, mainly depends 
on the composition of the interacting bodies and on the orbit. According to Newtonian 
Gravity Rr is approximately 

^(f)M^F^ < 23 > 

where p is the star average density. As usual, we define a penetration factor as the ratio 

P = ^. (24) 

where R p is the distance of closest approach between the two bodies. 

In order to establish what is the critical value of (3 below which the star can safely 
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approach the black hoel without being destroyed by the tidal interaction, we proceed as 
follows. We monitor the behaviour of the star principal axes as a function of time for in- 
creasing values of (3. As an example, in figure 1 we plot the axes for two parabolic orbits 
corresponding to f3 — 1 (upper panel) and f3 — 1.1 (lower panel log-scale). In the first case 
the axes oscillate around average values that approximately correspond to a conformal Rie- 
maniann stationary configuration (Diener et al. 1995). Conversely, for (3 = 1.1 we see one 
axis indefinitely growing, and this signals a destructive interaction. In table 1 we give the 
values of p cr u for circular, elliptic and parabolic orbits for the system we consider. (3 crit is 
the limiting value of (3 above which we observe the indefinite growing of one axis. 

In what follows we shall consider only values of f3 < (3 crit . In all cases we shall assume that 
the orbit lays on the equatorial plane ((x, y)-plane) and we shall compute the gravitational 
signal emerging in the z-direction in the TT-gauge. fo+ and will refer to the part of 
the signal due to the tidal interaction and to the consequent stellar deformation, computed 
using the quadrupole moment (21), whereas h+ b and h°£ b will be the orbital contribution, 
computed using eq (18). 



3.1 Circular orbits 

When the orbit is circular, the orbital contribution to the gravitational signal is 

h°;\t) = ih°; b (t) = A cos(2 ^ orbt) , (25) 

1 /2 

where u orb = (^jr) is the Keplerian angular velocity, D being the separation between the 
two bodies, r is the radial distance from the source, and 

A = (26) 

where M = M* + M B h is the total mass and \x = M * ^ BH is the reduced mass of the system. 
Thus, radiation is emitted in a single spectral line at a frequency v orh = i.e. 

h°;\u) = - IVorb) - S(u + 2v orb ] . (27) 

Due to gravitational emission the orbit decays; thus the signal we compute using eqs. (25) 
may not be correct, unless the changes induced on the orbit by the energy lost in GW occur 
on a timescale much longer than the orbital period, so that we can neglect radiation reaction 
effects. To check whether this is the case for the orbits we consider, we have computed the 
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time s 



Figure 1. The white dwarf principal axes are plotted versus time for two parabolic orbits with (3 = 1 and (3 = 1.1. The 
indefinite growth of one axis in the lower picture signals a destructive interaction. 



adiabatic timescale r 

E rb 



T ~ 



J GW 



(28) 



where E^b = — §^d* * s the system orbital energy, and L GW = dE £ t w is the GW-luminosity 
due to the time variation of the orbital quadrupole moment 
32 G 4 fi 2 M 3 



5 c 5 D 5 

The results of our calculations on the tidal interaction are summarized in figure 2 and 3 and 
in table 2. 
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Table 2. Orbital parameters and GW signals amplitudes are tabulated for different values of the penetration factor (3 (column 
1) for circular orbits. The orbital radius R, period T and frequency u or i, are given in columns 1 to r43; the adiabatic timcscale 
t is in column 5; in the last two column we give the amplitude of the spectral line of the orbital signal, A/ An, and that of the 
f-modc peak appearing in the tidal signal; they are both multiplied by the source distance, r. 



P R p 


T 


L'orb 




T 


(A/47r) x r 




10 4 (km) 


(s) 


(Hz) 




(s) 


10~ 4 (km) 


10~ 4 (km) 


0.15 14.6 


336.4 


3.0 x 10" 


-3 


4.8 x 10 11 


0.46 


1.7 


0.3 8.0 


118.9 


8.5 x 10" 


-3 


3.1 x 10 10 


0.87 


2.7 


0.4 6.0 


77.2 


1.3 x 10" 


-2 


1.0 x 10 10 


1.4 


4.2 


0.5 4.8 


55.3 


1.8 x 10" 


-2 


4.0 x 10 9 


1.5 


7.0 


0.6 4.0 


42.0 


2.4 x 10" 


-2 


1.9 x 10 9 


1.7 


9.2 


0.7 3.1 


33.4 


3.0 x 10" 


-2 


1.0 x 10 9 


2.5 


15 



In the upper panel of figure 2 we plot the Fourier transform of the tidal signal h+ (v) for 
(3 = 0.4. The plot shows a number of very sharp peaks. The first is at v — 2v or \ )) and it is a 
signature of the orbital motion on the tidal deformation of the star. The highest peak is at a 
frequency Vf = 0.146 Hz, which is the frequency of the fundamental mode of the non radial 
oscillations of the star in the unperturbed configuration (i.e. when it has a spherical form). 
This peak is surrounded by equally spaced peaks, at frequencies v = Vf±2nv or b, n = 1,2, 
due to the coupling of the orbital motion to the deformation. In addition, we also see smaller 
peaks at the frequencies of the first pressure mode, v Pl = 0.316, and further harmonics 
v — v pi ± 2nu or b. This picture clearly shows that the non radial modes of oscillations of the 
star are excited in the tidal interaction and their signature appears in the gravitational wave 
spectrum. The plot for h d ^ [y) is entirely similar and will not be shown. 

In the lower panel of figure 2 we show h+ {v) for increasing values of (3 (i.e. for decreasing 
orbital radii); we see that, as f3 increases and the orbit shrinks, the 2z/ or {,- peak shifts toward 
higher frequencies, but the peaks corresponding to the excitation of the f- and p-modes, of 
course, do not move. Moreover, due to che change of z/ or b the spacing between the harmonics 
changes as well. 

Finally, in figure 3 we show the signal for the case (3 = 0.7 which is close to tidal disruption 
(see table 1). In this case the contribution of different harmonics is higher, giving rise to a 
richer structure in the frequency spectrum. 

In table 2, for several values of (3 we show the orbital radius R, period T and frequency 
v or b (column 1 to 4); in column 5 we give the adiabatic timescale r. These data show that, 
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Circular Orbits 
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x 1e-05 
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0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 

v(Hz) 

Figure 2. The Fourier transform of the GW signal associated to the tidal interaction between a 1 Mq WD and a 10 Mq BH 
is plotted as a function of frequency for a circular orbit with /3 = 0.3 (up) and for different values of /3 (down). The main peak 
at vt = 0.146 Hz corresponds to the frequency of the fundamental mode of the non radial oscillations of the star, a smaller 
peak at v = 0.316 Hz indicates the excitation of the first p-mode. The remaining harmonics are due to the coupling between 
orbital and tidal motion (see text). 

for the orbits we consider, we are always in the adiabatic regime since r » T. In the last 
two column we give the amplitude of the spectral line of the orbital signal, A/Att, where 
A is given in eq. (26), and the amplitude of the f-mode peak appearing in the tidal signal 
h+. It is interesting to note that the f-mode occurs at frequencies higher than 2z/ or f,, and 
its amplitude is always larger than that of the obital signal. 

It is worth mentioning that the damping time of the fundamental mode is very large, of 
the order of tqw ~ 10 6 seconds for the stars we consider. 
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Circular Orbit 

v f 

- 2 v orb p=o.7 ; 

' 0.05 0.1 0.15 0.2 0.25 0.3 0.35 

v(Hz) 

Figure 3. Tidal signal for a circular orbit with j3 = 0.7. The coupling between orbital and tidal motion is stronger than in 
figure 2 (where /3 is smaller) and higher harmonics are excited. 

3.2 Elliptic orbits 

Whereas for a circular orbit radiation is emitted in a single spectral line at twice the orbital 
frequency, when the orbit is eccentric waves are emitted at frequencies multiple of z/ or b, 
and the number of equally spaced spectral lines increases with the eccentricity (Peters & 
Mathews 1963; Ferrari, D'Andrea & Berti 2000). In the left panel of figure 4 we compare 
the orbital (up) and the tidal (down) signal emitted on an orbit with f3 = 0.4 and e = 0.75. 

As in the circular case, the excitation of the f- and pi-modes is manifested by the sharp 
peaks at the corresponding frequencies, and it is interesting to see that the coupling between 
the orbital and tidal motion introduces a large number of harmonics. Note also that the 
maximum of the orbital signal is much larger than the f-mode peak in tidal signal. 

In the right panel of figure 4 we show how the structure of the tidal signal changes due to 
eccentricity, plotting h+ (v) for (3 = 0.4 and for two values of e, (upper panel): as expected, 
we see that lower eccentricities correspond to a smaller number of spectral lines due to the 
orbital-tidal coupling. In the same figure (lower panel) we compare the tidal signal emitted 
on orbits with the same eccentricty (e = 0.75) and different values of f3, showing how the 
mode excitation is sensitive to the penetration parameter. 

We find that for the orbital parameters tabulated in table 3, the maximum of the orbital 
signal is always much larger than the f-mode peak in tidal signal. 



19-04 
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0.045 



Elliptic Orbits 



Elliptic Orbits 




0.15 0.2 

v(Hz 

Figure 4. In the left panel we plot the Fourier transform of the orbital signal (up) and of the tidal signal (down) emitted by 
the WD-BH binary studied in this paper when the orbit is an ellipse with j3 = 0.4 and e = 0.75. As in the circular case, the 
main peak in h d ^ corresponds to the excitation o the f-mode, but its amplitude is much lower than that of the maximum in 
h°£ b . In the right panel we show the tidal signal for elliptic orbits with (3 assigned and varying e (up), and with e assigned and 
varying fi (down). 



Table 3. Orbital parameters of the elliptic orbits we consider, plotted for different /3 (colum 1): the eccentricity (column 2), 
pcriastron and apoastron distances (column 3 and 4), orbital period (column 5). 



13 


c 


R p 


Ra 


T 






(10 4 km) 


(10 5 km) 


(s) 


0.4 


0.25 


6.0 


1.0 


42 


0.3 


0.75 


8.1 


5.6 
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0.4 


0.75 


6.0 


4.2 


220 


0.5 


0.75 


4.8 


3.4 


157 


0.6 


0.95 


4.0 


15.7 
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3.3 Parabolic orbits 

While for circular and elliptic orbits the orbital GW-signal has a discrete structure in the 
frequency domain, parabolic orbits are associated to signals that are continuous functions 



the orbital and the tidal contributions are in a different frequency range; but now, unlike the 
circular and elliptic cases, the dominant contribution appears to be the tidal one, reaching a 
peak-amplitude of about ~ 1CT 2 at i/f = 0.146 Hz, about twice the amplitude of the orbital 
peak, which occurs at v — 1.6 • 10~ 2 Hz. 

Comparing the orbital signal of figure 5 with the orbital signal in figure 4 (upper panel, left) 
we see that the shape of the parabolic signal is basically the envelope of the elliptic one; 
indeed, the parabola can be seen as the limiting case of an ellipse with eccentricity equal 
to unity. It is interesting to compare the signal emitted on the parabolic orbit with that 
emitted on a very elongated elliptic orbit with the same (5: in the right panel of figure 5 
we show the orbital (up) and the tidal (down) signals emitted when the orbit is an ellipse 
with j3 = 0.6 and e = 0.95. We see that in this case the amplitude of the f-mode peak 
increases significantly with respect to the orbital one, though it does not become bigger. 
Thus it appears that, for a fixed /3, the oscillations modes of the star increasingly contribute 
to the total emitted radiation as the orbit elongation increases, and become dominant when 
the orbit is a parabola. 

In figure 6 we plot the tidal signal for parabolic orbits for different values of the penetra- 
tion factor. It should be mentioned that, as for circular and elliptic orbits, the signals emitted 
along the axis orthogonal to the orbital plane in the ' x '-polarization have a structure similar 
to that of the '+'- polarization, and for this reason we omit to show them. 

As a last example, we have computed the radiation emitted when the star penetrates more 
deeply into the black hole tidal radius, choosing values of j3 > 1, but always lower than the 
tidal disruption limit. In this case the star is highly deformed and assumes a 'cigar' shape; 
therefore, we cannot expect to see in the frequency spectrum the peaks corresponding to the 
excitation of the modes of the spherical star, as it was in previous cases. 
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Elliptic Orbit With High Eccentricity 




Figure 5. In the left panel we plot the Fourier transform of the orbital signal (up) and of the tidal signal (down) emitted when 
the WD moves on a parabolic orbit with f3 = 0.6. For comparison, in the right panel we show the orbital and the tidal signal 
when the orbit is elliptic, with the same f3 and very high eccentricity e = 0.95. 



However, the cigar also oscillates and his eingenfrequencies scale as the inverse of his length, 
being the mass fixed. Indeed the tidal signal presents a very interesting structure shown in 
figure 7. In the three panels on the left we show h d + } \v) for f3 — 1, 1.03, 1.05. We see that 
there is a dominant peak which corresponds to the excitation of the fundamental mode of the 
cigar; indeed, as expected, when (3 increases the tidal interaction is stronger, the star assumes 
a more elongated cigar shape and the frequency of the main peak decreases. The spectrum 
also shows several equally spaced spectral lines and their origin can be understood by looking 
at the right panel of figure 7, were we plot the principal axes of the star as functions of time. 
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Parabolic Orbits 
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Figure 6. The Fourier transform of the tidal signal emitted in parabolic orbits is plotted as a function of frequency for three 
values of /3. 

When the star approaches the black hole, one of the axes on the equatorial plane - let us 
name it 'a'- grows much more than the other two, 'b' and 'c'. In the figure the axes are 
normalized to their initial values, therefore at large distance (large negative time) they are 
all equal to 1 since the star is spherical. Analysing the axes behaviour, we can identify two 
characteristic frequencies: one, ui, which corresponds to the first oscillation, when the a-axis 
starts to grow, reaches the maximum elongation, and then decreases to the value about 
which it will oscillate in the new, approximately Riemaniann, stationary configuration; the 
second one, U2, corresponds to the oscillation of the a-axis about the new configuration. For 
the three values of (3 considered in figure 7 these frequencies are given in table 4. Referring, 
for example, to the signal h+ (v) for /3 = 1 in the left panel of figure 7, and indicating with 
vq the frequency corresponding to the fundamental mode of the 'cigar', we see that there 
are equally spaced spectral lines at v = v ± ni^, with n = 1,2,..., and smaller peaks (only 
the first on the left of the main peak is clearly visible) which correspond to v — i/ Q ± nvi . 
The spacing between the peaks corresponds exactely to the frequencies, given in table 4, at 
which the axes oscillate. 

3.4 Results for a different black hole mass 

It is interesting to see how the entire picture discussed above changes if we change the black 
hole mass. Be R p the periastron distance (coincident with the orbital radius in the circular 
case). The tidal tensor (5) scales as CV,- oc {Mbh/R^) = P 3 ', therefore we expect that, for an 
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Figure 7. The tidal signal emitted in a parabolic orbit is plotted for three values of /3 approaching the critical value (left 
panel). The green continuous line is the orbital signal, plotted for comparison. On the right panel we show the behaviour of 
the principal axes of the star as a function of time, when the star passes through the pcriastron. The relation between the axes 
behaviour and the emitted signal, and the origin of different peaks in h de f are explained in the text. 
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Table 4. In this table we give, for three values of j3 approaching the critical value and for parabolic orbits, three frequencies 
characteristic of the tidal deformation: uq, is the frequency of the main peak appearing in the tidal signals shown in figure 7 
and corresponds to the excitation of the fundamental mode of the deformed star; v\ , refers to the first oscillation of the a-axis 
as the star approaches the periastron, and corresponds to the oscillation of the a-axis about the new cigar-like configuration. 



P vo 


v\ 




(10- 2 Hz) 


(10- 2 Hz) 


(10-2 Hz) 
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Figure 8. The orbital signal emitted when the black hole mass is Mgjj = 10 s Mq is plotted for an eccentric orbit with j3 = 0.6 
and e = 0.75 (left) and for a parabolic orbit with /3 = 0.6 (right). These plots have to be compared to the upper panel of figure 
4 and 5, respectively, where the orbital signal is plotted for the same orbital parameters and a black hole mass Mgjj = XOMq. 

assigned (3, both the tidal signal h de f and the critical value of (3 needed for disruption are 
independent of the black hole mass. Indeed, this is what we find by numerical integration. 

Conversely, the orbital part of the signal, h orb , depends on the black hole mass. For 
instance, in the circular case from eq. (26) we see that for M* << Mbh the amplitude of the 
emitted spectral line is A oc Mbh, showing that the orbital signal increases with Mbh- A 
similar behaviour is exhibited in elliptic and parabolic orbit, as shown in figure 8, where we 
plot h+ b emitted when Mbh — 1O 3 M , for an elliptic orbit with (3 = 0.6 and e = 0.75 (left) 
and for a parabolic orbit with (3 = 0.6 (right). A comparison of these signals with those 
found for the same parameters when the black hole mass was Mbh = 1OM (upper panel, 
left, of figure 4 and 5, respectively) clearly shows how the orbital signal increases with Mbh- 
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4 CONCLUDING REMARKS 

In this paper we have computed the gravitational signal emitted when a white dwarf moves 
around a black hole on a closed or open orbit. We have computed both the orbital and 
the tidal contributions and compared the two, assuming that the star is close to the black 
hole, but in a region safe enough to prevent its tidal disruption. In all cases, the non radial 
oscillation modes of the star are excited (in our approach we do not include the dynamical 
behaviour of the black hole) to an extent which depends on how deep in the tidal radius the 
star penetrates and on the type of orbit. The orbital, h orh ', and the tidal, h def , contributions 
are emitted at different frequencies: smaller for h orh ', higher for h de} . 

The system we have analyzed in detail is composed of a 1 M Q white dwarf and a 10 M 
black hole. In this case we find that for circular orbits up to the critical one, the amplitude of 
the single spectral line in h orb is always smaller than that of the f-mode peak in h de *; however 
if the black hole mass is larger, the situation reverses, since h de f is nearly independent of 
M B h, while h orb is proportional to it. It is interesting to note that h def contains several 
peaks that shows the coupling between the orbital and tidal motion. 

When the orbit is elliptic, h orh shows several spectral lines emitted at multiples of the 
orbital frequency; h de f has again a main peak at the f-mode frequency, and in addition a 
number of peaks due to the orbital-tidal interaction, that contain a signature of the nature 
of the orbit. The amplitude of the f-mode peak of course increases with the penetration 
factor (3 and also, for a fixed f3, with the eccentricity; for the system we have considered 
this amplitude tends to that of the maximum of h orb for very elongated orbits (e > 0.95). 
Again changing the black hole mass to a higher value the orbital contribution increases with 
respect to the tidal one. 

For parabolic orbits the situation is very intersting, since the star can penetrate more 
deeply into the tidal radius, allowing for larger values of (5. As (5 tends to the critical value 
the star assumes a very elongated, cigar-like shape, deviating considerably by its initial 
spherical structure. Thus the main peak in h de f shifts toward lower frequencies, since the 
f-mode frequency scales as the inverse of the lenght of the principal axis. The harmonics 
which appear in the tidal signal corresepond to the oscillation frequencies of the principal 
axis of the 'cigar', and the larger is f3 the higher will be the number of excited harmonics. 
A comparison with the behaviour of the orbital signal (see figure 7) shows that if j3 < 1.03 
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the tidal signal is considerably higher than h orb if Mbh — 10 M , but again the situation 
reverses if the black hole mass is sufficiently higher. 

The results of our study show that the tidal signals emitted by a close encounter between 
a white dwarf and a black hole lay in a frequency region which is intermediate between 
the sensitivity region of LISA (10~ 4 — 1CT 1 Hz) and that of ground-based intereferometers 
VIRGO, LIGO, GEO, TAMA, sensitive above ~ 10 Hz. Detectors that fill this gap have 
recently been proposed, like the Big Bang Observatory, thought as follow on mission to 
LISA (Bender et al. 2005), and DECIGO, proposed by a Japaneese group (Seto, Kawamura 
& Nakamura 2001); they should both be extremely sensitive in the decihertz region 10 _1 — 1 
Hz, and would be the appropriate instruments to detect the tidal signals we have studied 
and to shed light on the dynamics of dense stellar clusters where these processes are more 
likely to occur. 
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